Photon bremsstrahlung and diffusive broadening of a hard jet 
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The photon bremsstrahlung rate from a quark jet produced in deep-inelastic scattering (DIS) 
off a large nucleus is studied in the collinear limit. The leading medium-enhanced higher twist 
corrections which describe the multiple scattering of the jet in the nucleus are re-summed to all 
orders of twist. The propagation of the jet in the absence of further radiative energy loss is shown 
to be governed by a transverse momentum diffusion equation. We compute the final photon spectrum 
in the limit of soft photons, taking into account the leading and next-to-leading terms in the photon 
momentum fraction y. In this limit, the photon spectrum in a physical gauge is shown to arise 
from two interfering sources: one where the initial hard scattering produces an off-shell quark which 
immediately radiates the photon and then undergoes subsequent soft re-scattering; alternatively 
the quark is produced on-shell and propagates through the medium until it is driven off-shell by 
re-scattering and radiates the photon. Our result has a simple formal structure as a product of 
the photon splitting function, the quark transverse momentum distribution coming from a diffusion 
equation and a dimensionless factor which encodes the effect of the interferences encountered by the 
propagating quark over the length of the medium. The destructive nature of such interferences in 
the small y limit are responsible for the origin of the Landau-Pomeranchuck-Migdal (LPM) effect. 
Along the way we also discuss possible implications for quark jets in hot nuclear matter. 



PACS numbers: 12.38.Mh, ll.10.Wx, 25.75.Dw 



I. INTRODUCTION 



The modification of high transverse momentum (high 
Pt) jets P, [3, [|[ as they pass through dense mat- 
ter has achieved center stage in the experimental pro- 
gram at the Relativistic Heavy-Ion Collider (RHIC) at 
Brookhaven National Laboratory (BNL) [J] and is poised 
for a more comprehensive exploration at the upcoming 
Large Hadron Collider (LHC). The central goal is an 
quantitative exploration of certain properties of the dense 
matter produced, based on the modification encountered 
by a jet while passing through the excited medium. The 
benchmark for such observables is provided by the mod- 
ification encountered by a hard jet in cold nuclear mat- 
ter The basic partonic processes which lead to the 
modification of the jet in cold or excited media are identi- 
cal and lead to a formulation in terms of the same corre- 
lation functions. Experimental observations of the mod- 
ification in cold and excited matter provide numerical 
values for these correlation functions and as such allow 
comparisons between these two environments. 

The most noticeable modification observed is the loss 
of forward momentum and energy of the hard partons 
believed to result from a combination of elastic [f| 0] 
and radiative [I]] mechanisms. The experimentally ob- 
served effect is the depletion of high momentum hadrons 
produced in the fragmentation of these partons after 
they escape the medium [f| HJ. For the light partons, 
the energy loss by induced gluon radiation is believed 
to be the dominant mechanism. Calculations based on 
this are by now considerably developed. There cur- 
rently exist four different schemes of radiative energy 



loss: the Arnold-Moore- Yaffe (AMY) [§], the Armesto- 
Salgado- Wiedemann (ASW) Gyulassy-Levai-Vitev 
(GLV) [ll| and the Higher-Twist (HT) []j approach, 
each incorporating a slightly different approximation 
scheme in the calculation of the radiative gluon spec- 
trum (There have also been attempts to calculate both 
the broadening of a jet and its energy loss in a dipole ap- 
proach [U[). The different schemes also utilize different 
definitions of the medium modified fragmentation func- 
tion and in most cases, apparently different models of the 
medium. 

The AMY and ASW approaches, based on the previous 
work of Baier et. al. (usually referred to as BDMPS) 0], 
are cast in the thick medium approximation, in the sense 
that the hard parton as well as the radiated gluon scatters 
multiple times in the medium and an infinite series of 
such scatterings has to be re-summed [HI, EH- Whereas, 
the GLV and HT approach have traditionally been cast in 
the thin medium approximation, where the hard parton 
and radiated gluon scatter a finite number of times in 
the medium, with the effect of each subsequent scattering 
being incorporated order by order in opacity or twist [lj| • 

Some time ago, it was pointed out that the higher twist 
formalism may also be extended to the thick medium 
limit [T3|. In a recent effort, the formulation of jet mod- 
ification due to multiple scattering within a re-summed 
higher twist approach was begun [l8|. In that paper, 
contributions from all twists, which are enhanced by the 
length of the medium, were re-summed to calculate the 
transverse broadening experienced by a hard parton as 
it traversed a dense medium without radiative emission. 
The resulting all-twist expression assumed the physically 
transparent form of a two dimensional transverse momen- 
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turn diffusion equation. The diffusion tensor was shown 
to be related, up to overall constant factors, to the well 
known jet transport coefficient q Q . It is the object of the 
current manuscript to extend this formulation of all-twist 
jet modification. The specific problem being addressed 
here is the spectrum of real photons radiated from a par- 
tonic jet multiply scattering off a dense medium. While 
the photon is assumed to escape the medium without in- 
teraction, very unlike the case of a radiated gluon, the 
problem does encode many of the characteristic features 
of gluon radiation such as the Landau-Pomeranchuck- 
Migdal (LPM) effect which suppresses very collinear ra- 
diation . 

Besides its role as an intermediate step to the prob- 
lem of energy loss by gluon radiation, hard photon 
bremsstrahlung from a high p? jet is a problem of suffi- 
cient phenomcnological relevance in it own right. Tradi- 
tionally, thermal photons and dileptons have been con- 
sidered as ideal probes of the plasma as they escape the 
plasma after their formation without any further interac- 
tion, directly conveying information regarding conditions 
at their origin [19| . Jet medium interactions, where a 
hard jet undergoes Compton scattering off a constituent 
of the medium leading to direct photon [2(| and dilep- 
ton 21( production have also been proposed as deep 



probes of the matter produced: these measure the par- 
tonic density of the medium at the scale of the hard par- 
ton. Recently, photon hadron correlations [22], [H, [24j 
where a hard photon is produced in the course of the mul- 
tiple forward scattering of a hard parton have emerged as 
an electromagnetic window to the energy loss transport 
parameter q. There now exist a series of measurements of 
photon and high px hadron correlations at RHIC. These 
are meant to be used in tandem with single inclusive high 
Pt hadron suppression [25| and hadron-hadron correla- 
tions [2(| to characterize the jet transport parameters of 
the medium. The results of the current manuscript are 
directly relevant to all near-side measurements of such 
correlations and represent an extra source of jet broad- 
ening in away side correlations between a trigger photon 
and a leading hadron. 

The results of this manuscript may also be easily 
ported into the calculation of hard photon production 
from a hot partonic medium. Such calculations, nowa- 
days, are mostly carried out within the finite tempera- 
ture formalism of AMY 0] by assuming that the medium 
exists at asymptotically high temperature and invoking 
the Hard Thermal Loop (HTL) approximation [27j (for a 
complete discussion see Ref . |28( and references therein) . 
Our reformulation of this problem in the HT approach 
achieves several goals. The AMY scheme invokes the as- 
sumption that the medium is a thermalized plasma of 
quarks and gluons at asymptotically high temperature 
T . The HT scheme makes no such approximation re- 
garding the medium; the only approximation made is 
that the medium have a short distance color correlation 
length. This is a property of a multitude of media includ- 
ing the HTL plasma. Hence, the results derived in the 



current article should be applicable to both cases where 
the HTL approximation is applicable and where it is not 
(for instance, in the confined environment of cold nuclear 
matter). This conjecture may be easily tested by calcu- 
lating the diffusion tensor in an hot partonic medium 
and promoting the hard parton from a single probe to an 
ensemble with a thermal distribution. In such a compar- 
ison, contributions that occur due to the initial parton 
being off-shell due to its formation in a violent collision 
should be neglected. 

Such comparisons and the associated phenomenology 
will be deferred to a future publication. Here, the focus 
will lie on a derivation of the all-twist re-summed ex- 
pression for photon bremsstrahlung from a hard parton 
traversing a dense medium. While the remaining sections 
specifically treat photon production from a hard jet pro- 
duced in the DIS off a large nucleus, it should be pointed 
out that this is merely an instrument of convenience and 
the results may be readily generalized to photon produc- 
tion off a hard jet in a different environment. To allow 
for such extensions in the future, our derivations are car- 
ried out in a factorized form, almost independent of the 
initial source of production of the hard quark. 

The paper is organized as follows: in Sect. II the sim- 
plest calculation of photon bremsstrahlung at leading 
twist is carried out. This is solely performed to famil- 
iarize the reader with the notation, frame and particular 
choice of gauge used in this manuscript. In Sect. Ill, 
the problem of photon radiation from a hard parton that 
scatters a finite number of times in the medium is ana- 
lyzed and general expressions are derived. In Sect. IV, 
sums over the various photon production points in the 
medium are carried out. The soft photon approximation 
is invoked and the twist expansion in cold nuclear matter 
is carried out in terms of a transverse momentum gradi- 
ent expansion. In Sect. V, the resulting expressions for 
soft photon production from a hard parton that scatters 
n times in the medium are re-summed to obtain the all- 
twist photon spectrum. Concluding discussions appear 
in Sect. VI. 



II. LEADING TWIST AND COLLINEAR 
PHOTON RADIATION 



In this section, the calculation of the photon produc- 
tion rate from a hard parton, produced in the DIS off 
a large nucleus, which subsequently undergoes multiple 
scattering in the confined nuclear environment is begun. 
The photon, not having scattered in the medium after 
its production, escapes the medium unscathed, revealing 
the nature of the matter at the location of its produc- 
tion. The interference between the photon produced at 
different locations in the hard parton's diffusive history 
lead to a modulation of the factorized rate of photon 
production. In this vein, we focus on an all twist eval- 
uation of the differential photon production rate from a 
hard jet. Consider the semi-inclusive process of DIS off 
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a nucleus in the Breit frame where one jet with a trans- 
verse momentum l q± and a bremsstrahlung photon with 
transverse momentum l± are produced, 



Ej jr , da 



dW^ 



d 3 L 2 d 2 l q± d 2 l ± dy 2ttsQ 4 ^ d 2 l q± d 2 l ± dy' 



(4) 



L(L 1 )+A(p)^C(L 2 ) + J(l q± )+ 1 (l ± ) + X. (1) 

In the above equation, L\ and L 2 represent the momen- 
tum of the incoming and outgoing lcptons. The incom- 
ing nucleus of atomic mass A is endowed with a mo- 
mentum Ap. In the final state, all high momentum 
hadrons [hi, h 2 , ...) with momenta pi,p 2 , . . . are detected 
and their momenta summed to obtain the jet momentum 
and X denotes that the process is semi-inclusive. 

Throughout this paper, the light-cone component no- 
tation for four vectors (p = [p + tP~,Pa_]) will be used, 
where, 



P 



P 



P ~P 



(2) 



The kinematics is defined in the Breit frame as sketched 
in Fig. [TJ In such a frame, the incoming virtual photon 
7* and the nucleus have momentum four vectors g, Pa 
given as, 



q = L 2 -Li = 



2q- 



• <1 



,0,0 



P A =A\p + , 0,0,0]. 



In this frame, the Bjorken variable is defined as xb = 
Q 2 /2p + q~ . The radiated photon has a transverse mo- 
mentum of lj_ and carries a fraction y of the forward 
momentum q~ of the quark originating in the hard scat- 
tering, i.e., 



y = — ■ 
q 



(3) 




Jdq) 



FIG. 1: The Lorentz frame chosen for the process where 
a nucleon in a large nucleus is struck by a hard space-like 
photon leading to the production of an outgoing parton and 
a radiated photon. 

The double differential cross section of the semi inclu- 
sive process with a jet with transverse momentum l qj ^ 
and a final state photon with transverse momentum l± 
may be expressed as 



where s = (p + L\) 2 is the total invariant mass of the 
lepton nucleon system. The reader may have already 
surmised the form of the leptonic tensor as, 



L 



(5) 



In the notation used in this paper, \A;p) represents the 
spin averaged initial state of an incoming nucleus with A 
nucleons with a momentum p per nucleon. The general 
final hadronic or partonic state is defined as \X). As a 
result, the semi-inclusive hadronic tensor in the nuclear 
state | A) may be defined as 



W A^ = Y^{^f8\q + P A -p x ) 



X 



x (A ]P \J^mx){X\J v mA;p) 



= 21m 



d 4 y e lqyo (A;p\J»(y )r(0)\A;p) 



(G) 



In the equation above, the sum (%2x) runs over all possi- 
ble hadronic states and J M = Qqipqj^ipq is the hadronic 
electromagnetic current, where, Q q is the charge of a 
quark of flavor q in units of the positron charge e. It is 
understood that the quark operators are written in the 
interaction picture, and two factors of the electromag- 
netic coupling constant a em have already been extracted 
and included in Eq. The leptonic tensor will not 

be discussed further. The focus in the remaining shall 
lie exclusively on the hadronic tensor. This tensor will 
be expanded order by order in a partonic basis with one 
photon in the final state and leading twist and maxi- 
mally length enhanced higher twist contributions will be 
isolated. 




FIG. 2: 



The Lowest order and leading twist contribution to 



The leading twist contribution is obtained by expand- 
ing the products of currents at leading order in a s and 
next-leading-order in a em to account for the radiated 
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photon. This contribution may be expressed diagram- 
matically as in Fig. [21 This represents the process where 
a hard quark, produced from one nuclcon in a deep- 
inelastic scattering on a nucleus, proceeds to radiate a 
hard photon and then exits the nucleus without further 
interaction. Other diagrams are suppressed in our choice 
of light cone gauge. In the following, we analyze this 
contribution in some detail. Indeed, there is no new in- 
formation presented in the current section and the dis- 
cussion of such contributions is now well established [2!| . 
The approximations made in this section will form the 
basis for the analysis of photon production at all twist. 
The semi-inclusive hadronic tensor may be expressed as, 



(7) 



AC* / d 4 y (p|V'(2/o)7 M 00 7XO)|p) 



C£ / d 4 y Tr[l-^^}F(y )O m (y ). 



In the above equation, Cp expresses the probability to 
find a nucleon state with momentum p inside a nucleus 
with A nuclcons. 

In the collinear limit, the incoming parton is assumed 
to be endowed with very high forward momentum (po + = 
xqp + ,Pq — > 0) with negligible transverse momentum 
Po± <C Pq ■ Within the kinematics chosen, the incom- 
ing virtual photon also has no transverse momentum. As 
a result, the produced final state parton also has a van- 
ishingly small transverse momentum (i.e., with a distri- 
bution S 2 (p±)). In this limit, the leading spin projection 
of the pieces which represent the initial state and final 
state may be taken. The factors, 







on 



Tr 



j_ 

2 

rf 4 / 



(: 



(10) 



x Tr 



(2^) 4 

7~ 



4 4 , d l q d po d p 
Z (2^) 4 (2^) 4 (2^) 4 



2 (p + q) 2 - ie. 



(Po + q) 2 + 1£_ 

x e n-yo e -i(po+q)-(yo-z) e —tt-(z—z') e —u q -(z-z') 
x e- l( - p '° +q > z 'e 2 q . 

Integrating over z and z 1 yields the two four-dimensional 
(^-functions: S 4 (po + q — I — l q ) and <5 4 (po — Po)- 

The erudite reader will have noticed that we have ig- 
nored various projections such as those which arise from 
the (-L)-components of the 7 matrices, e.g., 



C^Tr 



2 1 2 ' 



(11) 



This approximation may be justified in the high energy, 
collinear limit l\jy « Q 2 where such contributions are 
suppressed compared to those of Eq. (JT]). In so doing, 
the focus of the remainder of this article has been lim- 
ited to projections where the incoming virtual photon 
is transverse. As the primary impetus of the current 
manuscript is a description of the effect of final state 
scattering on the out-going quark leading to final state 
photon brcmsstrahlung, we defer the discussion of the 
slightly different process, where the initial virtual pho- 
ton is longitudinal, to a separate effort. 

The on-shell 5- function over I is used to set l + = 
l 2 L /2l^. The other on-shell <5-function, instills the con- 
dition, 



7° +7 3 



7 =7 



(8) 



are used to obtain the spin projections along the leading 
momenta of the outgoing state and the incoming state. 
The coefficients of these projections are the two func- 
tions, 



s(i 2 q ) = s[( Po + q-iy 



S[-Q 2 + 2p+{q- - r) - 2q+T - 2q-l+] 



1 



2p+q 
5[x - x B 



x (l - y)- x B {l - y) - 



/1 



2p+q y 



xl\ 



2p+q-(l-y) 



112) 



where, the collinear condition that po ± — > has been 
used to simplify the final equation. The new momentum 
fraction Xl has been introduced: 



F(y ) = A(p$(y o )±-ip(0)\p) 



(9) 



and (in a notation where the superscript on the operator 
00 implies that the quark undergoes no scattering in the 
initial or final state) 



XL 



2p+q y(l - y) p+T f 



(13) 



where y has already been defined as the momentum frac- 
tion of the radiated photon (l~ /q~) and Tf is the forma- 
tion time of the radiated photon. 
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The factor G a p(l) in Eq. (fT0|) represents the radiated 
photon's spin sum. In this effort, the light cone gauge 
(A~ = 0) will be assumed, i.e., 



dydl 2 . 



= C-*2nJ2Q 4 X(xB+x L )(-g^) 



n ■ I 



(14) 



2n I 



2 Pq-*n(y)- 



(17) 



where we have introduced the light cone vector n = 
[1,0,0,0] which yields I ■ n = l~. Note that with this 
choice of gauge, the largest component of the vector po- 
tentials from the initial states may still be regarded as 
the (+)-componcnts. 

Substituting the above simplifications in Eq. ()10l) . 
leads to the simplified form for the final state projection: 



VII) 



d 4 p dl~d 2 l 



x Tr 



(2tt) 4 (2tt) 3 2Z- q 
7" 



— e 2 e~ ipoVo 



2 2p+q (xo - xb - x Do - ie) 
x {tx7~([*o - x B ]p + - l + Wi{-9± a0 ) 

-7_L • <7 - 7 7x ■ I— rz— 



1- 



_2l+ 

7 7 + (<7 -l )l -rr 



7 9" 



x 2ir 



2p+<? (a; —x B — x D q + ie) 
6[x -x B - x L ] 



2p+q-(l-y) 



(15) 



The 5-function may be used to carry out the = XqP + 
integral; the absence of Pq and po ± from the integrands 
allows for these integrals to be carried out and constrain 
the locations yQ,yo ± to the origin. Further simplifica- 
tions may be carried out by noting that ^ anticommutes 
with 7_i_, while {7+, 7~} = 21 (i.e., twice the unit matrix 
in spinor space) and {7 =t ,7 ± } = 0. Replacing l~ = q~y, 



in Eq. (fT5|) . one obtains, 







(HI 



e 2 Q 2 A(q~) 2 



{2p+q-)Hp+q-(l-y)x\ 



2/ 2 2/+' 
2(x L p+-l+) + ^+2(l-y)fL_ 

1 v . 

Q 2 q a e m [dydl\ 2-2y 



-S(yt)s 2 (y 0± r q ( 



y 



(16) 



2tt J l\ y 

Rcintroduction of the final state projection O 00 in 
Eq. (|7|), produces the well known and physically clear 
formula for the differential semi-inclusive hadronic ten- 
sor with single photon emission in the final state, 



In the above equation, /^(ig) represent the parton dis- 
tribution function of a quark with flavor q and electric 
charge Q q in units of the electron charge e, in a nucleus 
with momentum Ap i.e., 



f?(xB+x L ) = A J ^L e -^B+^) P + 

x ltp\i'(yoh + m\p)- 



(18) 



In Eq. (17]), the factor P q ^ gy (y) = (2 - 2y + y 2 )/y is 
the quark-to-photon splitting function; it represents the 
probability that a quark will radiate a photon which will 
carry away a fraction y of its forward momentum. The 



projection = g 

As the parton produced immediately after the hard 
scattering has a vanishingly small transverse momentum, 
the transverse momentum of final produced quark is sim- 
ply the negative of the photon's transverse momentum, 
-l±. As a result, the differential hadronic 



i.e 



tensor for the transverse momentum distribution of the 
final quark is given as 



dW£ u 



dydl\d 2 l q ^ dydl 



dwr s 2 (i ± + i q±) . 



(19) 



III. 



PHOTON RADIATION FROM MULTIPLE 
SCATTERING 



In this section, the calculation of the single photon 
production rate from diagrams which include the mul- 
tiple scattering of the quark on the gluon field within 
the various nucleons inside the nucleus will be carried 
out. The focus will be on a generic diagram of the type 
that will be encountered in the evaluation of the photon 
differential rate at all twist. In another language, this im- 
plies the inclusion of higher twist corrections to Eq. (|17p . 
Higher twist contributions are obtained from diagrams 
which include expectation values of products contain- 
ing more partonic operators in the medium [30l ] e.g., 
the gluon field strength operator product F +v '(y)F+ (0). 
While contributions from the expectation of such oper- 
ators are suppressed by powers of the hard scale Q 2 , a 
sub-class of these are enhanced in extended media (such 
as an atomic nucleus) due to the longitudinal extent that 
must be travelled by the struck quark. As in Ref. [is| |. 
only such length enhanced higher twist contributions will 
be isolated and eventually rc-summcd. Issues relating to 
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the generalized factorization [3ll.l32j| of such contributions 
will not be dealt with in this effort; the focus will be to 
derive the effect of multiple scattering of the hard parton 
in the nucleus on its ability to radiate a collinear pho- 
ton. To obtain higher twist contributions, higher orders 
need to be included. A diagram with 2n gluon insertions 
may contribute to twist m < In. In all calculations, the 
high energy and hence small g limit will be assumed, as a 
result all diagrams containing the four-gluon vertex will 
be ignored. Within this choice of kinematics and gauge 
there exist a class of diagrams which contain three-gluon 
vertices and receive similar length enhancements. These 
diagrams involve radiated gluons in the final state and 
contribute to the energy loss of the parent parton. Such 
length enhanced contributions will be considered sepa- 



rately in an upcoming publication. 

The aim of this effort is to isolate the higher twist con- 
tributions at a given order in coupling and twist that 
carry the largest multiple of length L ~ A 1 / 3 [33|. In the 
language of power counting, one looks at the combination 
a™L n (usually n < m) and focuses on diagrams with the 
maximum n. The multiples of L are obtained by insisting 
on conditions that lead to the largest number of propa- 
gators going close to their on-shell conditions. This is 
easily achieved by the hard parton having the maximum 
number of space like exchanges with the medium. The 
photon is radiated on-shell and final state scatterings of 
the photon with the rest of the medium which involve 
additional powers of a em are ignored. 
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FIG. 3: An order m + n contribution to the single photon production rate. There are n out going gluons attached to the 
quark line on the right hand side of the cut and m incoming gluon lines on the left hand side of the cut. The photon attaches 
between the p th and (p + l) th location on the right and between the 5 th and (q + 1)* location on the left. This contributes to 
twist k < m + n. 



The generic contributions being considered, thus have 
the form of Fig. [3] A hard virtual photon strikes a hard 
quark in the nucleus with momentum p' (p in the com- 
plex conjugate) at location y' = (j/o m the complex con- 
jugate) and sends it back through the nuclear medium. 
At this stage, the quark has momentum q[ in the am- 
plitude and qi in the complex conjugate. In the course 
of its propagation, the hard parton scatters of the gluon 
field within the medium at locations, y[ with < I < m 
(at locations j/j with < i < n in the complex conjugate) 
wherein the hard parton picks up momenta p[ (pi in the 
complex conjugate) changing its momentum to q' l+1 (9i+i 



in the complex conjugate). The photon is radiated be- 
tween locations q and q+ 1 in the amplitude and between 
p and p + I in the complex conjugate. The quark line 
immediately after the photon radiation has momentum 
q' q+ 2 m the amplitude and q p +2 hi the complex conjugate. 
Following the notation of the previous section, the pho- 
ton has momentum / and the final cut quark propagator 
has momentum l q . The ellipses in between gluon lines 
in Fig. [3] are meant to indicate an arbitrary number of 
insertions. The various momentum equalities outlined in 
this paragraph may be now be expressed succinctly as, 
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q i+ i = q + ^Pt V[0<i<p] 
3=0 
fc 

9fe+i = ? + E^ V [0<fc< g ] 
z=o 

gi+i = q + /JPi - ^ V [p + 1 < i < n] 



(20) 



5=0 
k 



Vk+i = 1 + J2p'><- 1 V [g + 1 < < m] 



z=o 



The general hadronic tensor for such a contribution 
may be written as 



p,q J a 

A f ii 



(21) 



= AC* d^M^r^A 0»™ ^ 7 "V(0)|pi>. 

In the above equation, we have lumped all operators and 
propagators except those of the first struck quark in the 
amplitude and complex conjugate within the multiple op- 
erator product Op™. We have also suggestively factorized 
out one nucleon state denoted as \pi), which is the nu- 
cleon containing the struck quark. The remaining nu- 
cleons are contained in the state |^4) [34[ which is acted 
upon by the multiple operator product 0™ m . The coef- 
ficient Cp ± is meant to indicate the correlation between 
the first nucleon and the rest of the nucleus. We will 
defer issues related to such a factorization of the nuclear 
state to Sect. IV where the collinear expansion of the 
multiple operator product will be carried out and the 
nuclear state decomposed into an ensemble of nucleon 
states. The hadronic tensor, may be further decomposed 
in terms of the leading spin projections of the two oper- 
ator products, 



1 P-v 



■ct 



d%^\l-r\l v \F^)%^ -,(22) 



where, F(yo) is defined as in Eq. ^ with \p) replaced 
with \pi) and 0™'™ is the component of the leading spin 
projection of the final state operators, 



0! 



Tr 







p.q 



A 



(23) 



In the remainder of this section, we will focus exclusively 
on the component Op'™. 

The expression for 0™'™, may now be written down 
by applying the Feynman rules to the diagram of Fig. [3] 
The large rectangular blob is meant to indicate the nu- 
clear state. The lines connecting the quark propagators 
to the blob are not propagators themselves; these are sim- 
ply quark and gluon operator insertions, they introduce 
a certain momentum, spin and color into the diagram. 
Note that the quark operator insertions that couple with 
the initial incoming virtual photon have already been ex- 
tracted from the diagram and placed in Eq. (|2"Tj) . Note, 
there is no particular ordering of the gluon lines. As these 
are not propagators there is no meaning associated with 
crossed gluon lines. The entire set of n + m vertex inser- 
tions (with the gluon vector potentials contracted with 
the nucleus) may then be connected by quark propaga- 
tors in (n + m)\ ways (this overall combinatorial factor 
is removed by the (n + m) ! which appears in the denom- 
inator from the perturbation expansion). All the Feyn- 
man propagators arc written in position space as Fourier 
transforms of their momentum space expressions, e.g., 
for the case of a propagator in the complex conjugate 
amplitude between the i th and i + 1 th insertion, where 
i < p, we obtain, 



Tty(»W(tfc+i)] 



d 4 qi+i i 4i+ie 



-iqi+i-(vi-Vi+i) 



(2tt)< 
d 4 Pl *( 



-£j =0 ii) 



( 2 ^) 4 (9+e; =0 p,) 2 -^ 

x e -«(«+Ei=oPf )•(»*-»*+!). (24) 



Using the remaining relations between the various mo- 
menta as outlined in Eq. (|2"D|) the expression for the lead- 
ing spin component of final state operator product may 
be expressed as, 
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nrn _ TT TT d 4 y d 4 y' < ^ V% ~ 1 ^ V 'k-^ dH gjg 



P9 



i=l fe=i 



X Tl' 



1 P 

In 



(2tt) 4 (2tt) 4 (2tt) 4 (2tt) 



(25) 



i=0 



(E^o A)+ ^ / 
[(E;-= Pi) + # - « J ' t= P { [(Ej=o») + « - fl a - *e 



-7 



x 2^5(Z 2 )2^)G a/3 (Z) /, 



n r 



(Ef=o 



fc^-i I ' KEJUtf) +9-i] 2 +ieJ ' fe X 4 I [(EtoPi) + # + « J W + ^ + ie 

n— 1 m— 1 

n 1 



i=0 



X A 



k—m 



A 



The expression for 0™ Tl in Eq. (j25|) may be decomposed 
into a formal convolution of three terms, 

O™ 1 = J d 4 ld 4 l q T)y®y'®pT)p' 

x T(p > p / ,g > l > / g )r(p,p / > y > y , )M(y > y') (26) 



where, T denotes the trace in Eq. ([25]). which is solely 
a function of the momenta, M(y,y') = (A\ ■ ■ ■ \A) is 
the pure position dependent multi-operator matrix ele- 
ment and r(p, p', y, y') is the phase factor which contains 
both positions and momenta. The bold face quantities 
Pi y 5 p'j y' represent an array of momenta and positions, 



[Po, 



,Pn-l\] P 



[Pa, 



iPm-l]> 



[2/0, ■ ..,Z/nJ 



y' = [2/i,. 



The range of locations in the complex conjugate ampli- 
tude start from yg, which is the location of the initial 
struck quark, whereas, the location of the same parton 
in the amplitude is at the origin. Hence y' Q = is not 
explicitly written in the equations presented. The inte- 
gration measures S)p,p' and y' denote a product of 
integrals over the different four vectors contained in the 
arrays above. 

In an effort to simplify the writing of the matrix el- 
ement and the momentum dependent part, a certain 
aspect of collincar dynamics has been already been in- 
troduced. As calculations will be carried out in the 
light cone gauge (A~ = 0) in the Brcit frame at very 
high energy, the dominant components of the vector 
potential are the forward or (+)-components [33|, i.e., 
A a ~ g~" A + . As a result, the corresponding 7 matrices 
have solely a (— )-component. 



r 



To simplify the phase factor, an n th momentum (p n ) 
may be introduced, with 



1 = J d 4 Pn S 4 (l q -J2Pk-q + l\ 



(27) 



This leads to a considerable simplification of the phase 
factor as 



exp 



n m— 1 

J2 ■ w + w'j ■ y'j 

i=0 j=0 



i=0 



3=0 



(28) 



Note the complete absence of the initial hard photon mo- 
mentum q, or the final state momenta I, l q from the phase 
factor. The integrations over the lightcone components 
of the cut line l q have been performed using two of the 



four S- functions introduced in Eq. (|2 7p ■ The remaining 
two components of the transverse momentum integration 
(d 2 l q± ) arc constrained by thc(5-function, 



S 2 (C q± +l x -Kx) = S 2 f q± +l ± -^pi), (29) 



where, K± is a representative of the sum of the transverse 
momenta brought in by the n gluon insertions. 

The approximations stemming from collinear dynamics 
may now be instituted to further simplify the momentum 
dependent coefficient T in Eq. (|2"6"|) . The calculation is 
carried out in the Breit frame at very high energy. As 
a result, all momentum lines that originate in the target 
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are dominated by the large (+)-components of their mo- 
mentum, followed by their transverse coordinates, i.e., 



Pi >Pi± >Pl 



(30) 



A set of new quantities may now be introduced by or- 
ganizing the various factors which appear in the denomi- 
nators of Eq. (|25|) . There remain the obvious longitudinal 
momentum fractions, 



In most cases the above condition will allow us to prac- 
tically drop all (— )-components of momentum from for 
the hadronic tensor and focus solely on the (+) and (_L)- 
componcnts. It should be pointed out that the (— )- 
components are only being dropped from locations where 
they appear in addition to the larger (+, _L)-components. 
These are not dropped from the phase factors in T . 

The simplification of the numerator of T begins by iso- 
lating the largest components of the momentum. In most 
cases, this essentially reduces to retaining the sole factor 
7 + <7~ in the numerator of each of the propagators. A 
complication arises from the presence of a radiated pho- 
ton in the final state. The structure of the sum over 
polarizations controls which terms are to be retained. It 
should be pointed out that unlike the sum over polariza- 
tions in covariant gauge, 



G~ 



while, G" = — <7? p as in the covariant gauge. Here and 
in what follows, the _L tensor notation is introduced, i.e., 
A a ^" indicates that the only non-zero component of the 
tensor A involve its transverse components. The Dirac 
trace in the numerator of Eq. (|25|) denoted by N(T) may 
be separated from its propagation structure involving the 
denominator of the propagators and the ^-functions de- 
noted as D(T). These may then be evaluated separately, 
as in the following. 

The approximation is begun by ignoring all appear- 
ances of the (— )-components of the momentum which 
originate in the nuclear state i.e., p^ ,p[ , from all terms 
in the purely momentum dependent component of the 
integrand T, These momenta now appear solely in the 
phase factors, which allows for the p~ and p'~ integra- 
tions to be done, resulting in the localization of the pro- 
cess on the negative light-cone, i.e., 



0; 



a/3 



01+ I 1 - 
=l~- G ^+=G+^ = p,(31) 



r =1111/ d P i d Pi e != ° 

i=0 i=o 

12 iPi" (.Vi -V'ti) 

x e 1=0 

n m~ 1 

= n^ + -^)n^ + -^)- ( 32 ) 



i=0 



1=0 



In the above equation, y' = and as may be noted from 
the definition of the location arrays is not being inte- 
grated over. As a result, this constrains all the negative 
light cone locations in the above equation to the origin. 



Q 2 = 2x B p + q ; pt = x t p + ; p'j = x'jP + . (33) 

As in Ref. [l8| , we introduce two sets of momentum frac- 
tions which are dependent on the transverse momenta 
imparted to the struck quark from the medium, 



x Dl 



b Dk 



2p+q~ 



2p+q- 



(34) 



(35) 



Two new sets of momentum fractions which depend 
on both the transverse momentum imparted from the 
medium as well as that carried out by the radiated pho- 
ton are also introduced: 



3^0 

2y(l - y)p+q- 



\l±\ 2 -2/E2ZW± 
i=o . 

2y(i - y)p + q~ 



(36) 



(37) 



In the above equation, i, the index of the unprimed mo- 
menta (both for the longitudinal and transverse compo- 
nents) runs from to n, whereas, k, the index of the 
primed momenta runs from to m — 1, i.e., one less than 
maximum. The other indices, j and I denote smaller 
ranges of sums, with < j < i and < I < k. The reader 
will note the obvious role played by the factors above in 
the limit that the (— )-components of all momenta which 
originate in the nucleus have been ignored and one ex- 
pands around the dominant (-l-)-components along with 
the vectors q and I which have large (— )-componcnts. In 
this limit one may factor out the large components of the 
momentum i.e., p + ,q~ and express the remaining com- 
ponents such as combinations and products of l± , p]_ and 

p'_i_ as fractions of the large product p q such as XLi, 
xoi and x' Dk . 

Invoking the notion of the collinear approximation, the 
various definitions in Eqs. (|33H37p may be used to simplify 
the integrations that need to be performed in Eq. (|25|) . 
the factors appearing in D(T) may be written in a sim- 
plified form along with the longitudinal part of the phase 
factor (EHi as, 
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d(t) = n 



i=0 



2p + q ( X] _ XDl ~ x b - it 

3=0 



n 

i—p 



2p+q-(l-y) £>i~f 
J=0 



SB - £Li - if 



Vi=o 



2p+q-(l-y) 



n 

k—in — 1 



2p+«T (1 - y) - - x B - x' Lk + ie 



x exp 



v.J=0 
m — 1 



n 

k—q 
m — 1 



i=0 



1=0 



^i=0 1=0 



(38) 



In the above equation, the integration over x n may be performed with the (5-function which denotes the cut quark 
line to obtain, 



n-l 



(39) 



i=0 i=0 

This may be used to rearrange the longitudinal part of the phase factor as, 



exp 



E x Di 



1-2/ 



V 



p + (y n -y'm) 



J 



n-l 



ra—1 



J~J e [-iXiP + {Vi -Vn)] e [ ix ' l P + (y'l -»'m)] ) (40) 



-ix i p^(y i —y 
i=0 1=0 



where, the second line in the equation above involves only momentum fractions and locations from the left-hand side 
of the cut, where as the last line involves momentum fractions and locations from the right hand side of the cut 
line. The integrations over the remaining longitudinal momentum fractions, may now be performed starting from the 
propagators adjacent to the cut and proceeding to the propagators adjacent to the photon vertices. 

The first such integration, involves the propagator from the 3 rd line of Eq. ([55)1 . Isolating the piece that depends on 
the fraction x n -\ yields the integral, which may be performed by closing the contour of x~_ 1 with a counterclockwise 
semi-circle in the upper half of the complex plain of, 



dx n -i 



-ix ra -i-p + (!/„_i-3/ n ) 



Itt "-2 

Xn-1 + [Xi- X B - X Lr 

i=0 

i-y 



w(y n -Vn-x) e 



+ IB+IL„-1- E (^i-T^j 



(41) 



The effect of performing the above integration is the incorporation of both the real and imaginary parts of the above 
propagator into the overall expression of the hadronic tensor. It has the physical effect of propagating the quark from 
Vn—i t° Vn - Similarly, the integration over the propagator to the immediate right of the cut line may be carried out 
by closing the contour of x' m _ 1 with a clockwise semi-circle in the lower complex plain, : 
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Incorporation of the results of the above two integrals into the longitudinal phase factors leads to the expression, 



exp 



x exp 



' x Dn - Sx Ln 
k 1-1/ 



P Vn + * 



— 5x 



i-y 



n )p + v' m 



i\x B + XLn-l + — \ ° ] P + V n -1 +i\X B + X L ' m ~l 



n—2 m— 2 

]J cx p [- ix iP + (Vi ~ Vn-i)] II cx p ix 'ip + (y'7 - V'm-l) 

I 



i=0 



(=0 



Em— 1 / \ 
1=0 X Dl „+,/- 



(43) 



where, the notion of overall momentum conservation was 
invoked to define the new variable p' j_ = 22i=oP± ~ 

m 1 ~*k 



--Dm 



i=0 



XDi 



771—1 



(44) 



In this equation, the quantities labeled Xm for < i < 
n— 1 are defined in Eq. (f3~4")) , and x' D j for < j < to— 1 are 

defined in Eq. (|35|) . The n th transverse fraction i'o„ is set 
by the (5-function arising from the cut line and is given as 
in Eq. (|39p . The extra momentum fractions SxLm $ x 'Lm 
arise form the fact that there is a small difference between 

XLn and XLn-l i- e -> 



— 2Zj 



•ELn L n — 1 



•Pi* 



2p+q-{l-y) 1-y 



(45) 



and similarly, 



-21 



2p+q-(l-y) 



(46) 



It should be pointed out that Sxl u or the primed mo- 
mentum fraction has no dependence on the momentum 
fraction of the radiated gluon y. Thus such contribution 
will be retained even in the small y approximation. 

The evaluation of the remaining longitudinal momen- 
tum fraction integrals is analogous up to the outgoing 
photon vertices. Therefore, the general result for the in- 
tegrations over the remaining momentum fractions in the 
phase factor may be carried out up to the integral involv- 
ing Xp+i and Xq+i. Around the photon insertion there 
are two propagators, either of which may be evaluated 
to obtain the conditions on x p and x' q . We focus on the 
x p integral, the case for x' q is completely analogous. The 
integral in question is 



2?r 



p-1 



p "T /-I \ Xi ~ 1—y I ~ 1—y ~ X B ~ x Lp 



i=0 



p-1 



X P + Yl ( x i — x Di) - x Dp - X B - ^ 



iO(y p+ i - y p ) 

X L P + JZ^J 2-^=0 x Dj 



p + (y P -Vp+i) 



(47) 



The origin of two separate contributions lies in the fact that the leading length enhancement arises when at most one 
of the propagators is off-shell. This necessarily has to be one of the propagators adjacent to the outgoing photon line. 
The remaining integrations over the the propagators between the photon insertion and the hard scattering vertex, are 
similar to the case of transverse broadening [1 81 ] . 

The reader will have noticed that the integrals above are over the momentum fractions x% and x[, whereas, in 
Eq. (f25|) . the integrals are over the momenta, pf and p ; . We have substituted the definitions of Eq. |33|) and 
hence focus on the momentum fractions instead. Thus each integral produces a dimensionful factor p + . There being 
n + to + 1 integrals over the (+)-components of the momentum produces an overall factor of (p+)"+ m + x . This factor 
will not be explicitly written out in the remaining simplifications of the different parts of 0™'™ that follow. They 
will be reintroduced in the end of the section, when the different pieces are recombined to write down the simplified 
structure of 0™-™. 
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Invoking the above simplifications, the factor D(T), i.e., the part of 0™'™ which depends solely on the denominators 
of the propagators of the hard parton and the longitudinal phase factors in Eq. (|38[) takes the form, 



D(T) 



%r7 > Vn-1 > ■ ■ ■ > j/Q ) T] 

(2p+q-) n+1 (l - y) n ~P 



TT e -ix Di p+i/ t - J| e - i ( xn \-y'' , )p + Vi 



i=p+l 



x Lp + J2^=Q x Dj 



-i[xL p + T-^ T, x Di )p + y p -i{xL v + T-^ £aiDi )p + V p+ i 



%'m > y m-i > • ■ ■ > v'o ) 

(2p+q-) m+1 (l - y) m ~i 



l-\ 1 



k=0 



, lx 'D k P + v' k 



n - 

k=q+l 



M i-,, P' V k 



C Lq + l-j/ Ya=0 



2p+«-(l-y) 

In the above equation, the expression 9(y~ > y^-i > • • • ) is meant to indicate a product of n ^-functions: i.e., 



(48) 



%n > Vn-l >•■■) = II 



(49) 



i=0 



These ^-functions organize the longitudinal locations of the scatterings, starting from the hard incoming virtual 
photon vertex and moving towards the cut line (the same is true for the primed locations in the complex conjugate 
amplitude) . 

The numerator of the momentum dependent part of Eq. (|25p . which is denoted as N(T) may be decomposed, 
as in the previous section, into four parts [i.e., N(T) = ^2 i=1 N(T)i\ depending on the components of the photon 
polarization tensor chosen. These are denoted in the same order as in Eq. ([7]): N(T)i denotes the term arising 
from the projection G ±J ~. Using this combination of the components a, (3, the leading part of the numerator may be 
simplified as (in the remaining equations in this section, we will ignore the vector symbol on the transverse momenta 

l±iP±,p' j_; they should however, always be understood to be two dimensional vectors), 



iV(T)i=Tr 



/ v -\ 



\ II f 7" \ 7+T + 1 X XX \ ^« { ^(T - n + 7 X E^l - l± 



3=0 



9+1 



[ 7"7 + (T - n l-l + l-{-gf) 7"7 + (<T - n 



, i=p+l 



x 7 



\k—m—l 



v3=0 




r 



(50) 



NiT), = (2g-) 



— \n-\-m-\-l 



(i-y) 



n-Hm— p— q— 1 



In the above equation, the transverse components of commutation rule {7 ± ,7~ L } = 0, the entire numerator 
the propagators are only retained for those propagators term in the equation above may be readily simplified as, 
in the immediate vicinity of the photon insertion points, 
identified by the 7^,7^ in the above equation. Only 
the transverse components in these locations give leading 
contributions. The transverse components in the remain- 
ing propagators yield suppressed contributions compared 
to the longitudinal components of the momenta, q~ and 
q~ — l~ . Using the well known relations between the 
various 7 matrices, such as 7 ± 7 ± = 0, {7 ± ,7 =F } = 21 
(where 1 is the unit matrix in spinor space) and the anti- 



x 2 



p 



{i + d-y) 2 } E 



J=o 




E^+E/ 

.3=0 



1=0 



(51) 
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Using a similar procedure as above, the terms originating 
from the projection G +1 ~ and are written in combi- 
nation as (we refrain from presenting the full derivation 
of these terms as it is rather straightforward), 



N(T) 2 +N(T) 3 = {2q-) 



H n.\n+m-p-q 
- \n+m+l \ ± t) J 



x 2 



2ll-l±- EpI + E/x (2-|/) 



,i=0 



(52) 



Note that these factors appear solely due to the use of the 
light-cone gauge. In a covariant gauge such as Feynman 
gauge there are no photon polarization tensor compo- 
nents which connect the + and _!_ components. The last 
contribution originates from the projection G ++ , 



9/2 

N(T) 4 = (2q-) m + n+1 (l - y) m+n -P-i +1 2—± 
The sum of these four terms has the rather simple and physical form, 



(53) 



N(T) = (2q-) 



t (1 — 



2P 7 (y) 



\ix\ a -v\ix- E^l + E^ 



w=0 fc=0 



J=0 



(54) 



where, -P 7 (y) represents the photon splitting function. 

The full contribution to the term 0™>™ of Eq. (|23|) . post integration over all longitudinal fractions, may now be 
reconstituted by combining the various parts as stated in Eq. (|26j) : 



/• ; ,r- Mi " rti m ~ 1 i2/j 



2P 1 {y)e- lXBP v ° 

y(2p+q-) 2 (l-y) 2 

p 

n 

t=0 
9 

n 



x U e 

i=0 
9 

x Ue 

fc=0 



I Epj 



m 2 -w Efi+E^ 

t \i=0 k=o / J 



\l=0 



i=p+l 



ix' Dk P + y' k p-ip'x-v'x 



n 



1 II e (y'k -y'k-i) 



k—m 



x Lp + 1 _ y Yl P j=0 x Dj 
1 



X 'hq + 1-y Ya=0 X 'dI 



■^, + A;E4l JP »'« i^i g + r=? ; E ]p v'q+i 



X ( 



3 " i+ "Tr 



3=n 



A;p 



(55) 



In order to simplify the discussion of the remaining sections, the soft photon radiation approximation, i.e., y — > 
will be introduced. Only the leading behavior in y will be retained. This leads to a considerable simplification of the 
various expressions and also represents a well studied limit. Two of the energy loss formalisms, those of Rcfs. (Io| 
and [ll| are cast in this limit. In the limit of the soft radiation approximation, all factors of 1 — y in the denominators 
may be replaced with 1 + y in the numerators and the leading and next-to-leading contributions in y retained. As 
a result of this approximation, only the contributions from N(T)4 and N(T)2 + N(T)s need to be retained as they 
contains the leading and sub- leading negative power of y and are thus dominant in the limit y -C 1. In this limit, 
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|;j 2 (l + y)-y£22±-pl 



3=0 



bll 2 



2p+q y 



> yxDi 



E2pi -pi 

J'=0 



(56) 



As a result, retaining solely the leading and next-to-leading behavior in y we obtain 0™^", in the soft radiation limit 



■J2i jj- J2; ™ m " j2„t J2„/ 

o^" = e 2 / 1 i tt d \\ n d »rAi n n n ^^m^+zw?^ 



P,9 



X 



2P 1 (y)e- lXBP y o 



(2tt) 3 2Z- (2tt) 2 



i=0 v ' j=Q y ' 



y(2p+q-) 2 (l-y) 2 



\i=0 fc=0 



n < 

i=p+l 



-i(x Di -Sx Li )(l+y)p+y i Jp\-y\ 



q m 1 

ix 'n k p + y'k p-ip'l-y'l TT p i ( s r>k-*« ; £k)(i+i/)p + i''k P -Vrvl 



X 

fe=0 



1 r e -"ipP + y P _ e -«i.pP + yp + i 

XL-r, L 



fe=g+l 

1 



x ( A;p 



g m+n Tr 



1 



e ixl L q p + v' q _ e ixl L P p + y' q +i~ 



2—1 j = Tl' 



A;p 



(57) 



r 



The expression derived above is completely general, in 
the sense that no assumption regarding the nature of 
the nuclear state has been made. In the next section, 
a re-summation over the different locations of where the 
photon is radiated will be carried out. Following which, a 
factorization of the above hadronic tensor into a part that 
is solely dependent on hard momenta and a part depen- 
dent on soft momenta will be carried out and simplifying 
assumptions regarding the nuclear state made. This will 
then be followed by the re-summation over number of 
scatterings. 



IV. SUM OVER PHOTON PRODUCTION 
POINTS, FACTORIZATION AND GRADIENT 
EXPANSION 



out over the various quantities p, q and n, m. A few com- 
ments are in order: in analogy with the case of transverse 
momentum broadening of a hard quark propagating in a 
nuclear medium [l8| , cases where n ^ m represent unitar- 
ity corrections to terms where the produced quark scat- 
ters min[?i, m] times. The cases where p = q represent 
squares of the amplitude where the photon is produced 
in the p th location. The cases where p ^ q represent 
interference terms which lead to the well known Landau- 
Pomcranchuck-Migdal (LPM) suppression of the photon 
production rate. 

In this section, the sum over the various locations of the 
radiated photon in the amplitude and the complex con- 
jugate amplitude (i.e., p, q ) will be carried out. Invoking 
the small y approximation, the l± dependent momentum 
fractions may be expanded as, 



In the preceding section, a general expression for the 
single photon bremsstrahlung from a multiply scatter- 
ing hard quark parton was derived. The expression in 
Eq. (j57) represents the case where the quark scatters 
n times in the amplitude and m times in the complex 
conjugate. The photon is produced after the p th scatter- 
ing in the amplitude and after the q th scattering in the 
complex conjugate. To obtain the full differential cross 
section to produce a photon with momentum fraction y 
and transverse momentum l±, sums have to be carried 



xl p = x L (l + y) ~y^Jx Li 



i=0 



2p+q~y 2p+q- 



(58) 



As most 



To simplify notation, we set Ax l L = Ylj—p 5x Lj ■ 
of the factors in 0™^ n [Eq. (|57|) ] are independent of p, q, 
this contribution may be written compactly as, 
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y(2p+q-) 2 (l-y) 



2 



i=0 j=Q 

where, the terms that depend on the location where the photon is produced are denoted separately. The product of 
these factors is given, in the small y limit as, 



i - y + 



Ax p 



XL 



i-y 



Ax'? 



XL 



i[(x L (l+y)-Ax p L ]p + y p +i£ Sx Li p + y i -i[x L (l+y)-Ax p L }p + y p+1 + Sx Li p + y i 
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In both the denominators and in the exponent, the leading contribution is from the term xl ~ l/y- The remaining 
terms Ax 1 ^, Sxli an d x % D are all independent of, hence subleading in, y. Therefore, the factor of y in products such 
as 6xLi(l + y) may be dropped. 

The sums to be carried out are those over < p < n — 1 and < q < m — 1 restricted to the terms ,f p g q . The 
reader will note that these two sums are independent of each other and as a result may be carried out rather easily. 
The sum over p is 
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.(60) 



Note that the 2 nd phase factor in the p th term is identical to the I s * phase factor in the (p + l) th term. As a result, 
many terms cancel in the sum if the small correction Ax p L in the coefficient of the phases is neglected. The phases 
which never cancel are the I s * phase of the term p = and the 2 nd phase of the term p = n — 1. The result of the 
sum over p, up to next-to-leading terms in y, is 
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(61) 



In the limit that the photon is extremely soft, one may retain solely the leading contribution in y, i.e., only the first 
two terms in the brackets above. This corresponds to the "deep-LPM" limit: The regime where the radiated photon 
cannot resolve the various scatterings and treats the entire nucleus (or the entire process of n scatterings) coherently 
as one single scattering event. The next-to-leading term in y which is the third term in the equation above represents 
the first correction to this limit where the photon just begins to resolve the different scatterings of the hard quark 
in the medium. The terms obtained from the sum Y^=o m ^ nc complex conjugate amplitude yield a near identical 
expression with n replaced by m and the phase factors complex conjugated. 

We may now add the terms f n and g m i.e., the cases where the photon is produced after all the n or m scatterings. 
Such contributions necessarily have an on-shell, cut quark line on one side of the photon vertex. As a result, these 
terms do not contain the difference of two phase factors as in Eq. (|Tf|) and only have the first factor which occurs 
when the line after the photon emission is taken on shell (in this case taken by cutting the line). For the sum over p 
this yields, 
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We will also ignore the next-to-leading terms in y that 
occur in the phases as these will not play a major role in 
what follows. With these simplifications, we obtain 
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(63) 



Up to this point, the collinear approximation has been 
used to simplify the expressions for the hadronic tensor, 
without the introduction of factorization. The separa- 
tion of the hadronic tensor into a hard short distance 
piece and a soft long distance contribution may now be 
accomplished. All factors in Eq. (|57[) which contain the 
hard scales p + ,q~ constitute the hard part. All factors 
that depend solely on the soft _L momenta and distances 
along with the matrix clement constitute the long dis- 
tance element. All phase factors which contain a factor 
x l D p + (as well as its primed counterparts) or xlP + as part 
of their arguments belong in the hard part. The purely 
transverse phase factors such as exp[ip^_ • y\} (which do 
not contain any factor of l±) belong in the soft part along 
with the matrix elements. Thus, we may decompose the 
general contribution to the hadronic tensor as, 



£ 



0;;T= I ®yVp±H( P + ,q- 7 p ± ,l ± ,y)S(p ±7 y).(64) 



In the above equation, y and p± arc representative of the 
entire set of distances and transverse momentum that 
appear in Eq. (|5"T|) . 

The soft part which contains the matrix elements of 
the gluon vector potentials in the nuclear state may be 
simplified first. The simplifications arise from approxi- 
mations made regarding the nature of the nuclear state 
\A; p). In this manuscript, the nucleus is approximated as 
a weakly interacting homogeneous gas of nucleons. Such 
an approximation is only sensible at very high energy, 
where, due to time dilation, the nucleons appear to travel 
in straight lines almost independent of each other over the 
interval of the interaction of the hard probe. In a sense, 



all forms of correlators between nucleons (spin, momen- 
tum, etc.) are assumed to be rather suppressed. As a 
result, the expectation of the n + n' + 2 operators in the 
nuclear state may be decomposed as 



n+n 
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(65) 



where, the factor C px represents the probability to find 
a nucleon in the vicinity of the location y, which is a 
number of order unity (it is the probability to find one 
of A nucleons distributed in a volume of size cA within 
a nucleon size sphere centered at y). The remaining co- 
efficients Ct\ represent the weak position correlations 
between different nucleons. The overall factor of A arises 
from the determination of the origin (the location in 
the equation above) in the nucleus, which may be sit- 
uated on any of the A nucleons. Solely for the current 
discussion, we reintroduce the quark operators 4>(y~ , y±) 
and ^(O) in the above equation. 

It is clear from the above decomposition that the 
largest contribution arises from the term where the ex- 
pectation of each partonic operator is evaluated in sepa- 
rate nucleon states as the yi integrations may be carried 
out over the nuclear volume. As a nucleon is a color 
singlet, any combination of quark or gluon field strength 
insertions in a nucleon state must itself be restricted to 
a color singlet combination. As a result, the expectation 
of single partonic operators in nucleon states is vanish- 
ing. The first (and hence largest) non-zero contribution 
emanates from the terms where the quark operators in 
the singlet color combination are evaluated in a nucleon 
state and the n + n' gluons are divided into pairs of sin- 
glet combinations, with each singlet pair evaluated in a 
separate nucleon state. This requires that n + n' is even 
and may lead to a maximum overall factor of 



C 



^[(n+n'+2)/2] ; 



(66) 
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in the large A limit. It should be pointed out that large 
contributions may also arise, in principle, when n + n' 
is odd. In this case, the two quarks and a gluon are 
considered in the singlet combination with the remain- 
ing gluons evaluated in singlet pairs in the remaining 
nucleons. Here we institute the experimental observa- 
tion that (n)-parton observables are much smaller than 
(n — l)-parton observables. This is only true, once again, 
outside the saturation regime. In this effort, the focus 
remains exclusively outside this region, as a result we ig- 
nore all terms with more than two quarks or two gluons 
per nucleoli. 

Further simplifications arise in the evaluation of gluon 
pairs in a singlet combination in the nucleon states by 
carrying out the y± integrations. The basic object under 
consideration is (ignoring the longitudinal positions and 
color indices on the vector potentials) 



d 2 yid 2 y" ± (p\A+(y\)A+(y' J ± )\p) 



= (2n) 2 5 2 (pi -p'{) [ d 2 y^e-' lx '^ + ^- y '^ 



x ^ (p\A+ (y x /2)A+(-y ± /2) \p) , 



(67) 



where, y x is the transverse gap between the two gluon in- 
sertions and p± = (pj_ +pj_)/2. The physics of the above 
equation is essentially the transverse translation symme- 
try of the two gluon correlator in a very large nucleus. 
One will note that the two dimensional delta function 
over the transverse momenta has removed an integration 
over the transverse area of the nucleus thus reducing the 
overall A enhancement that may be obtained. This is 
then used to equate the transverse momenta emanating 
from the two gluon insertions in the amplitude and com- 
plex conjugate amplitude. This also simplifies the lon- 
gitudinal phase factors which now depends solely on the 
difference of the longitudinal positions of the two gluon 
insertions. In order to simplify further discussion, we 
consider the specific case where there are n gluon inser- 
tions in the amplitude and m = n + 2 An insertions in the 
complex conjugate (where An is a positive integer). Un- 
der the assumptions of restricting each nucleon state to 
be acted upon by only two gluon operators, we find that 
n + An nucleons are involved. The n path ordered gluon 
insertions in the amplitude are matched up by a series 
of path ordered insertions in the complex conjugate am- 
plitude on the same set of nucleons. The remaining An 
nucleons with 2 gluon insertions each may be distributed 
between the n nucleons. The equating of transverse mo- 
menta emanating from and being injected into these An 
nucleons insists that they impart vanishing net trans- 
verse momentum to the propagating quark. As such, 
the transverse momenta emanating from such nucleons 
does not contribute to the overall transverse broadening 
of the quark and so does not have any appearance in the 
overall transverse momentum delta function. Thus the 



overall transverse momentum ^-function only contains n 
different transverse momentum variables from p\ to p\ . 

One now invokes the collinear approximation in ex- 
panding the hard part as a Taylor expansion in trans- 
verse momenta around the origin p\ — ► 0. In the case of 
a cut with m = n + 2An there are 2n gluon insertions 
which produce transverse momentum broadening and as 
a result, as many derivatives, which involve n different 
transverse momenta. All terms of the form 
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if 



i p H \ P ±=oP± P± > 
i=i 1 dp ± dp ± 

where, p < n, yield gauge corrections for the contri- 
butions with 2p gluon insertions and as many deriva- 
tives [33|. The genuine 2n twist correction at this order 
has p = n. Expanding to this order, one obtains the 
generic term, 
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Hi 
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tr^pi>f H(P '^'^' } U=o 

x piVx^i (wr.fi/2Xj 0/7 > -vi/ 2 )- ( 68 ) 

Where, we have assumed the result of Eq. (f6"T|) and rein- 
troduced the color indices and longitudinal locations. 

Using integration by parts over the transverse distance 
y\ one may convert the product p\ A+.(y~ ,y 1 ± /2) — > 
^d"A^.(y~ ,y l ± /2) . In the extreme collinear limit, in 
the presence of a hard scale such that g is small, one 
may make the approximation, 



d±A+ ~ F a + 7 



(69) 



where, F a \_ represents the gluon field strength. Carrying 
this out consistently on the two gluon operator in the 
nucleon state of Eq. ([67]) and ignoring derivatives of the 
field strength \d a F +l3 ~ g(m 2 )g al3 j + — > 0] we obtain, 



rf 2 yxPl Q pl /3 e^-^(P^ + (^/2)A+(-y ± /2)|p) 



d 2 y ± e^^-(p\F +a (yj2)F + P(~yj2)\p) 
= f d 2 y^y-^(p\F+P{yj2)F+{-y^/2)\p){70) 



In the last line of the above equation we have averaged 
over the spins in the two field strength expection in the 
nucleon state with the constraint that the operator being 
evaluated in the nucleon be a spin singlet. The nucleon 
states are always assumed to be spin singlets or in spin 
averaged states. 

With the derivative expansion (at vanishing transverse 
momenta pj_ — ► 0) imposed on the hard part H, it no 
longer has any functional dependence on the transverse 
momenta. The integrations over the transverse momenta 
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may now be included completely into the soft part. As in 
the case of transverse broadening [l8[ , the action of the 
transverse momentum derivatives on the phase factors 
will not be considered. All such derivatives necessarily 
extract a factor of spatial separation yj (with < j < 
n, m), as 
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(71) 



These result in spatial moments of the two gluon matrix 
elements such as (p\F +p (y~ ,y J j_)y~ F p +(0)\p), which will 
be ignored in this effort. 



The action of the 2n derivatives is thus restricted to 
the overall transverse momentum conserving delta func- 
tion and the transverse momentum factors of Eq. (|63|) . 
The remnant transverse momentum dependence lies in 
the soft part which contains the matrix elements of the 
gluon field strengths and (5-functions which equate pairs 
of transverse momentum. We point out that, while we 
are considering cases where n = to — 2 An the final results 
for n = to + 2An are identical. The entire structure of 
the term with n + to gluon insertions and 2n transverse 
derivatives may be expressed as (we also approximate 
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In the above equation, n ~ 1 PAn represents the number of permutations of An pairs of consecutive gluon insertions 
in n — 1 locations. The barred transverse momenta p± = (p± +p' ± )/2, whereas Sp' l± = p' l± — represent the 

consecutive gluon pair insertions with two gluons per nucleon. The transverse momentum 5-functions in the soft 
part that correspond to the second type of insertions require that the momentum brought in by one such gluon is 
taken out immediately by the other. As a result, such insertions play no role in the double differential transverse 
momentum distribution being calculated and may be ignored. They however provide unitarity corrections to the 
total cross section. Alternatively, one may state that such insertions induce two exactly opposing currents in close 
proximity and thus produce a vanishing net transverse broadening and completely destructive radiative contributions. 
In the remainder we will consider only symmetric contributions, where n = to or An = 0. 

In the symmetric case, one may carry out p'\ integrations using the transverse momentum (5-functions in the soft 
part. As a result, there is a pair-wise equality between the transverse momentum brought in by the gluon insertions 
on the left and right hand side of the cuts. Using integration by parts one may also replace d/dp'\ — > —d/dp\. 
Setting n = ?n in Eq. (|72p . the action of the derivatives d/dpi± (where 1 < i < n) on the transverse momentum 
dependence in the hard part may be simplified as, 
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(73) 



where the transverse derivatives may act entirely on the transverse momentum 5-function or one or two derivatives 
may act on the terms within the curly brackets above and the remaining may act on the <5-function. The p\ derivatives 
acting on the transverse momentum <5-function may be replaced as 
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Using the above equation, the action of the 2n derivatives on the combination of the two dimensional transverse 
momentum dependent ^-function and the phase factor ^ p f p g q followed by an imposition of the limit of very small 
transverse momenta p l ± — > may be simply expressed as, 
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The action of the derivatives on the transverse 6- 
function or on the factors of p\ that appear in f p g q are 
the only non- vanishing contributions from Eq. (|57p . The 
longitudinal integrals, due to color confinement, yield the 
requirement that the longitudinal locations of the two 
gluons which act on the same nucleon state be in close 
proximity. The factor of Sy~ represents the small gap 
between the longitudinal positions between the gluon in- 
sertions in a single nucleon. One now tries to identify 
the most length enhanced term by isolating the maxi- 
mum number of unconstrained dy~ integrals. Note that, 
due to the assumption of short distance color correlation 
(ignoring color and spin indices), 



dy-dy' (p\F(y-)F(y' )\p) ~ dy~{FF)y-, (76) 
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This reduces the overall trace over color factors to 
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(77) 
(78) 
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The remaining n longitudinal position integrals for the 
gluon insertions may be simplified as 



where, (FF) is the gluon expectation at the mean lo- 
cation y~ ± y~ and y~ represents the color correlation 
length in the medium. In a nucleus, this is equivalent to 
the confining distance, whereas in a quark gluon plasma 
it would be related to the Debye length. 

Each such integral yields a factor of L~ ~ A 1 ' 3 from 
the unconstrained y~ integration. Equating the pairs of 
transverse momenta that appear in each two-gluon cor- 
relation, as well as using the relation between the lon- 
gitudinal momenta from the 9- functions in Eq. (|5T|) . re- 
quire that the largest length enhancement arises from the 
terms where the gluon correlations are built up in a mir- 
ror symmetric fashion, i.e., where the gluon insertion at 
y % is contracted with that at y' % , One may now let the 
transverse momenta p\ (for all i) in the hard part tend 
to zero and integrate over the remaining p\ integrals in 
the soft part, setting the corresponding transverse dis- 
tances between the gluon insertions in a nucleon to zero 
i.e., y\ -> 0. 

One may average colors of the quark and gluon field 
operators, 
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j UdyrOfo -id) = - { J \{dyr. (80) 



Invoking the above simplifications, the leading length en- 
hanced contribution at order 2n to the term O nn is ob- 
tained as, 
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In the above equation, the function D, represents the 
following expectation value in a nucleon state, 
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5 = tr^ |i?o+a(y " )Fa ^ (0)|p) ' (82) 

While no limits have been specified in the integration over 
longitudinal separation y~ , it should be understood that 
due to the assumption of a short distance color correla- 
tion length in the medium [as in Eq. I|76|) ]. the integral in 
the above equation only receives strength form the region 
with \y~\ < y c . The nuclear functions, £ ,± represent the 
following expectation value in the nuclear state, 

[ L) ~ 2N C J 2tt 

x c OF (p\F a+a {y; +y-)F a +( Vi ;)\p).(83) 

As in the case of the correlator D, the integration over 
the longitudinal separation y~ is limited by the color cor- 
relation length in the medium or the confining distance 
in a nucleus. The integration over the mean location of 
the two gluon field strength insertions i.e., y~ has its up- 
per bound constrained only by the size of the medium L 
(in the case of a nucleus L ~ A 1 / 3 ). The lower bound 
of y~ , in the case of E + , is given by y~: the location 
of the ip operator insertion in Eq. ([55)) . In the case of 
E~ , the lower bound is y~ = 0, given by the location of 
the ip operator in Eq. (fo5[) . It should be pointed out the 
two lower bounds are separated by at most one unit of 
the color correlation length (the size of a nucleon in the 
case of DIS on a large nucleus) and as such, this small 
distinction will be ignored in what follows. 

While the two nucleon states on either side in the above 
equation are denoted as carrying the same momentum p, 
this may not strictly be the case. Throughout this cal- 
culation, the hard scattering piece, which includes the 
scattering of the initial virtual photon with the incom- 
ing quark (including the nucleon which contained the 
quark) was factored out [see Eq. (|2"2"j) ]. In the case of 
the first term in Eq. (|8Tj) , which represents the case that 
the photon is radiated immediately after the hard scat- 
tering, such a factorization along with the deconvolution 
of the nucleus into nucleon states is well justified as there 
is no longitudinal momentum which is picked up in the 
multiple scattering which leads to the transverse broad- 
ening. In the case that the photon is radiated at a later 
scattering, the nucleon state which is struck by the hard 
virtual photon and the nucleon state which radiates the 
final outgoing photon may exchange longitudinal momen- 
tum differently between them in the amplitude and in 
the complex conjugate. In such a case, these two nu- 
cleon states represent the convolution of two generalized 
parton distribution (GPD) functions (37|, where the ex- 
cess momentum from one such distribution is balanced 
by the other. The momentum correlation between these 
two states, no doubt, depends on the distance y~ between 
them. The overall constant cof in Eq. (|83[) accounts for 
both the average momentum correlation in a large nu- 
cleus and corrections resulting from the approximation 



of replacing a generalized parton distribution by a reg- 
ular diagonal parton distribution. A phenomcnological 
discussion on cof and its value in large nuclei may be 
found in Refs. [38[. In the next section, a re-summation 
of all such contributions of arbitrary order n will be car- 
ried out. 



V. RE-SUMMATION AND SOFT PHOTON 
PRODUCTION FROM MULTIPLE SCATTERING 

In the preceding section, the length enhanced contri- 
butions of twist-2n to the soft photon double differen- 
tial rate were extracted. It turned out that the pho- 
ton rate at twist-2n (including the leading and next- 
to-leading contributions in y) depended on two different 
correlation functions D and E evaluated in the nuclear 
state, as well as on the multi-nucleon combinatorial fac- 
tor C^ o „ . In the case where terms such as these are 
not very small compared to the leading twist contribu- 
tion, all such terms need to be re-summed. 

In order to sum over all n, this last factor has to be 
simplified. In a sense, re-summation requires that this 
coefficient have a formally multiplicative structure. In 
the preceding section, a model of the nucleus as a weakly 
interacting homogeneous gas of nucleons was used. The 
formal expression of this assumption is hidden within 
the dimensionful parameter Cp p „. The precise eval- 
uation of such combinatorial coefficients is rather com- 
plicated, even for the case of next-to-leading twist [38| . 
From general dimensional arguments, in the case of non- 
interacting nucleons, this factor may be approximated 
as, 

C£L.^*C?(^p) n . ( 84 ) 

where, p is the nucleon density inside the nucleus and 
l/2p + originates in the normalization of the nucleon 
state. The remaining unknown coefficient Cp is now 
considered to be independent of the order n and is in- 
cluded in the leading twist hadronic tensor of Eq. (7}. 
This decomposition is also identical to that used in the 
case of transverse broadening. It should be pointed out 
that there may still be an extra normalization factor for 
each nucleon and this may be included with the factor 
of p/2p + . We will ignore any such normalization factors, 
as they are not essential to the re-summation. It should 
however be pointed out that the correlation between the 
nucleons off which the hard parton scatters and the origi- 
nal nucleon struck by the hard virtual photon is different 
from that between the original nucleon and that where 
the outgoing photon is emitted. As mentioned above, 
these constitute off-forward distributions where there is 
momentum shared between the two nucleon states. 
The decomposition of the combinatorial factor allows 
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for the definition of the new quantities, 

D = ^ C = W'/^ |F " + " F °»' (85) 

and similarly, 

E ± (x L ) = ^E ± . (86) 

Note that D is exactly the diffusion tensor for the trans- 
verse broadening experienced by a hard parton travers- 
ing the nucleus without radiation. This is directly pro- 
portional to the well known jet transport parameter 




new coefficient c p accounts for the weak correlation be- 
tween the nucleon struck by the hard virtual photon and 
the nucleon from which the outgoing photon is radiated. 
In principle, c p may depend on the shared momentum 
fraction xl- 

Both the terms in the above equation may be re- 
summcd with the observation that the sum, 

°° (dl~v? Y 

Hi q± ,L-) = £ a — -^ § 2 (i ± + r q± ), 

obeys the diffusion equation, 



q = 



8D. 



(87) 



The new quantities i? ± may be related to D as, 

E ± (x L ) = J d yp -coFe ±lXLP+ ^- y ° } D( yp ;). 

It should be pointed out that D is a local quantity evalu- 
ated at particular location y p and could very well change 
with time in a non-static medium, whereas the £?'s are 
intcrgrated over the space time path of the jet as it tra- 
verses the medium. 

In the following, we will only present the re-summation 
of the quantity 0™'", i.e., 



= E ° n ' n - 



(89) 
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The re-summed differential hadronic tensor and the dif- 
ferential cross section may be easily obtained by incor- 
poration of the above quantity into Eqs. (|4l7p . The final 
state operator matrix element for a parton undergoing 
1 < n < oo scatterings may be re-expressed as, 
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(90) 



In the above equation, it is understood that the second 
term in the square bracket does not receive contributions 
from the case of n = 0. The reader may immediately 
verify that for the case of n — 0, which corresponds to 
the case of no scattering, the first term by itself repro- 
duces the result of Eq. (fTo| for the case of y — ► 0. The 
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4{l q± ,L-) = DVi6{l q± ,L-). 



(91) 



with the initial condition discerned from Eq. (fT!)|) as, 

4>(L- = 0, l q± ) = 6 2 (l q± +T ± ). (92) 

As demonstrated in the case of transverse broadening, 
contributions such as those from Eq. (|72|) lead to a uni- 
tarization of the cross sections and, as a result, normalize 
the solutions of the diffusion equation above. The general 
normalized solution to Eq (f9Tj) is given as , 



<KL~,lx) 



1 



AnDE 



■ cxp 



l± +L 



4DL- 



(93) 



The reader will note, that (j>(L~,l±) reconverts back to 
the two dimensional delta function in the limit of L~ — > 
and the solution is unitary in the sense that 







dL- 



J d 2 l ± (t){l ± ,L-) ~ 0. 



Substitution of the above solution back into Eq. (1901) , 
yields the very simple expression for the final state re- 
summed operator matrix element 0, i.e., 



= 



2vr 



dfidy 



d 2 l q± P 1 (y)e- i[xB+XL)p+v » 



1 + yc p 



{E+(x L ) + E-(x L )}l'i + l ± -l q _ 



2DL- 



x 4>{L ,l q± ). 



(94) 



Substitution of the above final state operator matrix el- 
ement into Eq. (J7J), yields the full differential hadronic 
tensor and as a result, the differential cross section for 
photon production at all twist. 
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dW 



dydl\d 2 l q ^ 



2tt 



;2 



2tt 



e -i{ XB +x L )p+y Fgfo) 



1 + yc p 



{E+(x L )+E-{x L )} 



lj_ + l± ■ lq_ 

T 2 " 



2DL- 
<j>(L-,l g± ). (95) 



The above equation represents the main result of this 
article. The unintegrated expectation value of the two 
quark operator F(yo) is defined in Eq. ©. The differ- 
ential spectrum of photons radiated from a hard parton 
undergoing multiple scattering in the medium has been 
expressed as a factorized product of the photon split- 
ting function, the two dimensional transverse momentum 
distribution of the propagating parton and a multiplica- 
tive factor [factor in square brackets in Eq. (|95p ] which 
encodes the phases picked up by the particles as they 
traverse the medium. It is this factor which has to be 
evaluated for each path taken by the hard quark in the 
medium, weighted by the expectation of the two gluon 
field strength operator product along the path. The un- 
known coefficients, Cp and c p have to be determined by 
a phenomenological analysis. The decomposition of the 
above equation into quark and gluon structure functions 
in the nucleus and its associated phenomenological im- 
plications will be discussed in an upcoming publication. 



VI. CONCLUSIONS AND DISCUSSIONS 

The study of the modification of hard jets in dense 
matter is now approaching a rather sophisticated stage 
of its development. The wide variety of experimental 
measurements of the modification of jets and jet like cor- 
relations require the emergence of a single formalism ca- 
pable of describing the modification of hard jets over a 
wide range of energies and a variety of different media. In 
this manuscript, the higher-twist mechanism of jet mod- 
ification is extended in an effort to make it applicable to 
both the thin medium and thick medium limits. The lim- 
its of thick and thin refer to the number of scatterings 
that a hard parton will encounter on its path through 
the medium. While previous efforts [H, [l6[ have focused 
on the modification due to a few scatterings, the current 
manuscript furthers the development of the formalism 
required to describe the modification due to an infinite 
number of scatterings. 

The focus of the current manuscript has been to com- 
pute the spectrum of single photon brcmsstrahlung from 
a propagating parton undergoing multiple scattering in 
an extended medium. The entire process is cast in 
a Deep-Inelastic Scattering framework, where a virtual 



photon (with virtuality Q 2 ) strikes a hard parton in a 
nucleon, which is itself situated in a large nucleus. The 
parton, on its way through the nucleus, scatters multiple 
times and radiates a real photon. Energy loss due to real 
gluon emission was ignored. Each scattering of the hard 
parton is suppressed by powers of Q 2 , but is enhanced 
by the length of the nucleus ( cx A 1 / 3 ) over which the 
scattering may occur. For large A, an infinite class of 
such contributions need to be re-summed. 

The photon may be radiated from the vicinity of any 
of the scattering locations. The sum over such contribu- 
tions, leads to the destructive interference for very for- 
ward radiation known as the LPM effect. In the extreme 
case of an on-shell parton entering such a medium and 
radiating a very soft photon with momentum fraction 
y << 1, it is shown that the destructive interference is 
complete. This is referred to as the deep-LPM limit. In 
the current effort, the first two corrections in terms of 
initial virtuality and expansion in powers of y, to this 
limit, are computed. The first correction is due to the 
fact the the initial parton that enters the medium is not 
on-shell but may indeed be very virtual. As a result, it 
may radiate the photon even without the need for fur- 
ther re-scattering. The second contribution originates 
from not assuming the very soft limit y — > for the pho- 
ton, instead keeping the first leading corrections in y. 
This is the reason our results differ from the case of no 
re-scattering by a factor proportional to y. 

As in the case of transverse momentum broadening, 
it is demonstrated that the infinite series of power cor- 
rections may indeed be re-summed. The result has the 
simple and intuitive structure of a product of the photon 
splitting function, the two dimensional transverse mo- 
mentum distribution of the propagating parton and a 
multiplicative factor which encodes the phases picked up 
by the particles as they traverse the medium. As in the 
case of multiple scattering without radiation, a Gaussian 
profile in transverse momentum centered around the ori- 
gin was obtained for the propagating parton [l8[. How- 
ever, the Gaussian obtained in this case was found cen- 
tered around the momentum — Zj_, needed to balance the 
transverse momentum of the photon. 

The results of such a computation are crucial to the un- 
derstanding of both the near side and away side correla- 
tions between a trigger hadron and an associated photon 
or vice- verse. The phenomenological comparisons with 
experiment and other theoretical approaches, based on 
the expressions calculated above, will appear in a future 
publication. The derived results also bear considerable 
relevance secondary source of single hard photon 

production. 

An obvious next step is the evaluation of gluon radi- 
ation in the same process. Ultimately this might lead 
the way to a unified description of multiple scattering in 
the medium and radiative energy loss. While the photon 
does not scatter in the medium, its formation incorpo- 
rates many similar physics issues such as the destructive 
interference of the LPM effect. Such effects will reappear 
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in the future calculation of gluon radiation and will have 
to be dealt with in the case of energy loss of hard partons. 
However, the various scales that will be encountered in 
that problem as well as the approximation scheme used 
will be near identical as that in the current manuscript. 
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